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Abstract. Using several cosmological observations, i.e. the cosmic microwave background anisotropics 
(WMAP), the weak gravitational lensing (CFHTLS), the measurements of baryon acoustic oscilla- 
tions (SDSS-i-WiggleZ), the most recent observational Hubble parameter data, the Union2.1 compi- 
lation of type la supernovae, and the HST prior, we impose constraints on the sum of neutrino masses 
(Z f^v), the effective number of neutrino species (A^eflf) and dark energy equation of state (w), indi- 
vidually and collectively. We find that a tight upper limit on ^ niy can be extracted from the full data 
combination, if A^eff and w are fixed. However this upper bound is severely weakened if A^gff and w 
are allowed to vary. This result naturally raises questions on the robustness of previous strict upper 
bounds on ^ niy, ever reported in the literature. The best-fit values from our most generalized con- 
straint read 2 = 0.556+°-23i eV, A^eff - 3.839 ± 0.452, and w - -1.058 ± 0.088 at 68% confidence 
level, which shows a firm lower limit on total neutrino mass, favors an extra light degree of freedom, 
and supports the cosmological constant model. The cunent weak lensing data are already helpful in 
constraining cosmological model parameters for fixed w. The dataset of Hubble parameter gains nu- 
merous advantages over supernovae when w = -I, particularly its illuminating power in constraining 
A'^eff. As long as w is included as a free parameter, it is still the standardizable candles of type la 
supernovae that play the most dominant role in the parameter constraints. 
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1 Introduction 

One of the most intriguing mysteries in modem cosmology is the property of cosmic neutrinos. 
Directly from the oscillation experiments [1], it is known that 

AmL - lAm^J = (2.2!°;^) x lO'^eV^ 

AmL - Amj, = (8.1!";^) x lO'^eV^, (1.1) 

where Amj. = mj - rn^. is the squared mass difference between two neutrino eigenstates. This result 
naturally leads to three scenarios for the mass splitting of the standard three flavor neutrinos [2] : 

(I) Normal Hierarchy — mi ~ 0,m2 ~ A/Mgun, and nii, ~ Amatm; 

(II) Inverted Hierarchy — mi ~ m2 ~ Amatm > and m^ ~ 0; 

(III) Total Degeneracy — mi ~ mj ~ m^ » Aniatm- 

Although the absolute values of neutrino masses are beyond the reach of oscillation experiments, 
other ground-based experiments, e.g. measurements of tritium beta decay [3] and neutrinoless dou- 
ble beta decay [4], provide some hints on the overall mass scale. Nevertheless, the most com- 
pelling bounds on the total mass of neutrinos (^ my) come from the cosmos (for recent reviews, 
see Refs. [5, 6]). Firstly, the anisotropics of the cosmic microwave background (CMB) are sensitive 
to neutrino masses [7]. As massive neutrinos become non-relativistic, their contributions to the en- 
ergy constituent changes from radiation-like to matter-like. Provided that these neutrinos are massive 
enough, they will slow down to non-relativistic speed even before the recombination epoch. This 
slight modification in expansion history affects acoustic peaks. Small scale anisotropics are mean- 
while boosted, due to the enhancement of photon energy density fluctuations by the diminishing 
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gravitational potentials [8]. The latest observations of WMAP seven-year pin down the 95% upper 
limit of niv to 1.3 eV [9], while Planck and other next-term probes are presumably capable of push- 
ing forward the sensitivity to ^ my ~ 0.2 - 0.03 eV, [10-12]. Secondly, massive neutrinos affect the 
evolution of matter perturbation in a particular way: they erase the density contrast at wavelengths 
smaller than their characteristic free-streaming scale, resulting in a suppression of linear matter power 
spectrum as given by [13], 

APM ^ _3 

with Oy = 2 my/ ^93.8 eV^ being the neutrino matter density fraction [2]. Thus the constraint on 
the sum of neutrino masses greatly benefits from accurate measurements of matter power spectrum. 
A great deal of effort has been dedicated to reducing the upper limit of ^ my from a number of 
large scale structure projects, e.g. 2dFGRS, SDSS, WiggleZ, CFHTLS galaxy surveys [14-32]. The 
present record gives ^ wiy < 0.26 eV at 95% confidence level (CL), using the combination of WMAP 
seven-year data, SDSS DR8 LRG angular spectra together with the HST prior [33]. 

Another novel and independent tracer of matter clustering is weak gravitational lensing (WL), 
i.e. the bending of light through intervening inhomogeneous mass distribution, which as a result 
generates small distortion (at the 1% level) on the images of distant sources. Since the first detection 
of cosmic shear [34-37], WL has been repeatedly shown to be a powerful and precise approach to 
measure mass fluctuations and constrain cosmological models (see Ref. [38] for a thorough review). 
The greatest advantage of WL is that it only relies on the total matter content along the line of 
sight, and is therefore free from the problematic modeling of galaxy-to-matter bias [39]. Using 
Fisher formalism, several studies have forecasted the error budget of future WL surveys constraining 
neutrino masses [40^3]. However, there are few groups who have put actual constraints on ^ my 
from current observational results [44]. In particular, using up-to-date WL data of CFHTLS-T0003 
and other cosmological probes, Ref. [45] and Ref. [46] obtained the mutually consistent result, an 
95% upper bound of ^ my < 0.54 eV. 

Known as cosmic chronometers, the observational Hubble parameter data (OHD) have recently 
gained increasing attention in their potential of measuring the geometry and matter content of the uni- 
verse [47, 48]. The Hubble expansion rates at different redshifts can be obtained via the differential- 
age technique [49, 50]. Compared with the type la supernovae (SNIa) data, OHD have the virtue of 
being directly linked to the expansion history, whereas what SNIa measure is luminosity distance, 
related to H{z) through integration. Using the latest OHD data, Ref. [51] has derived competitive 
upper limits on the total neutrino mass, ^ my < 0.24 eV at 68% CL. 

Although there has been a vast scope of work on the constraints on ^ my previously, most of 
these analyses do not seek appropriate treatments of correlations associated with ^ my and other 
cosmological parameters. As argued by Ref. [52], the degeneracy between ^ my and the dark energy 
equation of state (EOS) parameter w fatally weakens the constraints on ^ my. Except the role of 
dark energy, the effective number of neutrino species, i.e. A^eff , remarkably lifts the bounds on ^ my, 
attributed to their correlation [53-55]. Standard primordial nucleosynthesis with weak interaction 
rate corrections yields A^gfT = 3.046 [56]. Here this non-integer part originates from the non-thermal 
feature in neutrino distribution funtion due to partial heating during e- annihilations, including both 
finite temperature QED corrections and flavor oscillation effects [57-59]. However, as proven more 
and more firmly by the series of WMAP papers, A^eff seems to be much larger than the standard 
value [9, 60, 61]. Similar conclusions are also reached by other small scale CMB observations, e.g. 
A^eff ^ 4.6 ± 0.8 for ACT, and A^eff ^ 3.86 ± 0.42 for SPT [62, 63]. The presence of extra light degrees 
of freedom, called "dark radiation", might be reasonably explained by some light sterile neutrinos, 
axions or majorons in thermal equilibrium after the QCD phase transition [64-66]. As a consequence. 
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most of the upper limits on X ^v, which are derived fixing A^gff - 3.046, seem too optimistic. We 
note that recently Ref. [67] uses several updated cosmological probes to deliver one of the most 
comprehensive investigations on the niy, Nf^s, w) constraints and their correlations. However 
Ref. [67] does not focus on comparing different sets of cosmological data, and furthermore it does 
not use one of the most promising techniques, i.e. WL. In this work, we would like to concentrate 
on the correlations between some key pai^ameters including my, N^e, w), and also investigate 
how their individual values and degeneracies affect the constraints as well as dark matter clustering. 
Particularly, we will study the role of WL, and compare the constraining ability of SNIa and OHD. 

The outline of our paper is as follows. Sect. 2 introduces the observational data in our analysis. 
Then we present the main constraint results in Sect. 3. In Sect. 4 we investigate in detail the impacts 
of the three key parameters (X my, Nes, w) on structure formation. Finally Sect. 5 presents our 
conclusion and discussion. 



2 Observables and Data 



In this section, we briefly introduce the cosmological probes selected for neutrino constraints. In the 
context of Friedmann-Robertson-Walker metric, it is simple to derive the coordinate distance r in 
terms of the radial distance as [68] 



X 

c 



sinh(M;^) if f2k>0, 
if - 0, 

^^x] if a < 0, 



(2.1) 
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with X and r also known as the comoving distance and the comoving angular diameter distance 
respectively [69]. Here the comoving distance at a certain redshift z' is given by 



X{z) = \ 
no Jo 



dz 
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where E{z.) denotes the expansion rate of the universe [70], i.e 
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with a = 1/(1 + z) standing for the scale factor. Into the above formulae, the property of cosmic 
neutrinos enters via 



- ^DM + fib + ^y> 

\4/3 



(2.4) 
(2.5) 



Besides Qy and A^eff> here Q.^, Odm, ^^b, ^r, ^y, ^k, corresponds to present-day density frac- 
tions of all the matter, cold dark matter, baryon, radiation energy, photon energy, curvature, and dark 
energy, respectively. Hq is the Hubble constant and c serves as the speed of light. 

Cosmic massive neutrinos directly influence the expansion history of the Universe. They also 
leave their signatures on structure formation and mass distribution, which can be investigated through 
large scale structure and WL data. To use these data, one needs to know the exact theoretical pre- 
diction of the matter power spectrum, i.e. P{k), defined as the Fourier transform of the two-point 
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correlation function of the matter density contrast 6{r) - I - h{r)/n{r). The linear perturbation theory 
provides reliable results of P{k) at large scales, usually referred to as the linear regime [7 1-73] . How- 
ever, as the wavenumber k rises up toward ~ 0.1 - 0.6 h/Mpc, the linear perturbation theory starts to 
break down, and various effects of non-linear clustering take place [74]. The scale-invariant suppres- 
sion on P{k) caused by massive neutrinos (Eq. 1.2) thus no longer works. Fortunately, people have 
developed numerous ways to explore those small scale phenomena, e.g. the one-loop perturbation 
theorem [75-77], which usually takes up the third order expansion of P{k) as a close approximation, 
the time renormalization group flow method [78], and the most widely quoted HALOFIT analytical 
formulae as the mapping algorithm calibrated by N-body simulations [79], 

A^(^) = ^ = + a2 , (2.6) 

where the total power spectrum A^(k) is split into a quasi-linear term Aq and a pure non-linear term 
A^. However, three problematic issues are found to aff'ect the neutrino damping effect reproduced by 
this method [80]: 

(I) HALOFIT overpredicts the suppression in non-linear regime, with maximum discrepancy 
occuring at ~ 1 h/Mpc; 

(II) HALOFIT cannot capture the redshift dependence of maximum suppresion's location in 
power spectrum; 

(III) at very small scales (k > 2 h/Mpc), HALOFIT underpredicts the non-linear power. 

Ref. [80] conducted a more extensive suite of N-body simulations particularly designed for 
quantifying precisely the suppression from massive neutrinos, and improved the original method into 
HALOFIT-v algorithm, 

A'\k) = {a'I + A'y • (1 + Qy), (2.7) 

with both quasi- and non-linear terms modified, and Qy accounts for additional non-linear growth of 
the neutrino component [cf. Ref. 80, appendix A]. For the sake of the accuracy of our constraints, we 
adopt HALOFIT-v as the correction to the original Pim{k) in presence of massive neutrinos, calculated 
numerically through the hierarchical Boltzmann equations. 

2.1 Weak Gravitational Lensing 

After the calculations of matter power spectra, under Limber's approximation, the convergence power 
spectrum of WL can be computed from the integration of matter power spectra along line of sight 
(see e.g. Ref. [81]), 

where I - k • rix) represents the multipole moment and «(y(z)) accounts for the source redshift 
distribution. In principle, xn should cover the entire Hubble radius, yet we can still safely replace 
it by ^(ziim)> i-C. the comoving distance out to the survey limited redshift (zum), where the source 
normalized number count (?i(z)) is sufficiently decaying. Since massive neutrinos damp matter power 
spectrum at small scales, the convergence power of WL is also reduced. 

For the actual survey, we choose the published CFHTLS-T0003 observations which have mea- 
sured about 2 X 10^ galaxies with /ab magnitudes between 21.5 and 24.5, imaged on an area of 57 
square degree (35 square degree effectively) [82]. CFHTLS-T0003 has provided the community with 
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the observational results of the shear correlation function ^e,b> the shear top-hat variance (lyP)gg> 

and the aperture-mass variance ^M^p^. Amongst these three types of WL measurements, (^MIJ^ has 
the virtue of less systematics prone due to its unambiguous local decomposition [83]. Theoretically, 
the aperture-mass variance is related to the convergence power spectrum via [84, 85] , 

K)(^^ = ^^l^PU^^^P'^ (2-9) 

with Jaix) being the Bessel function of the first kind. Due to the choice of WL data, we use the 
following parameterization for the source distribution, as recommended by Ref. [82], 

niz)oc-^, (2.10) 
z" + c 

r-ziim 

I n{z)dz = 1, 
Jo 

where a, b, and c are nuisance parameters ' and the survey limited redshift is assigned the value of 
Ziim = 7 (M. Kilbinger, private communication). 

2.2 Type la Supernovae 

The use of SNIa as standardizable candles provides a powerful probe to explore the properties of dark 
energy. SNIa observations furnish one of the metric distances, i.e. the luminosity distance, 

Di^iz) = r{z){l + Z) (2.11) 

where r(z) is the comoving angular diameter distance given by Eq. 2.1. We use the latest Union2.1 
complation of SNIa dataset reported by Ref. [86] in the redshift range 0.015 <z< 1.414. This SNIa 
catalog consists of 580 individual supernova events and is cautiously calibrated against numerous 
sources of systematic uncertainties. Compared with the former Union2 compilation [87], this updated 
dataset includes twenty-three new events at high redshifts (0.6 < z < 1-4) and thus are helpful in 
tightening the constraints on the early behavior of dark energy. 

Besides the Union2.1 compilation, there are a number of popular SNIa samples, e.g. ESSENCE 
[88], Constitution [89], SDSS-II [90], SNLS3 [91]. In deriving those datasets, different hght curve fit- 
ters are adopted. Generally speaking, three independent fitters are proposed and maintained amongst 
the cosmology community, i.e. SALT2 [92], MLCS2k2 [93], and SiFTO [94]. When trained on the 
same SNIa data, SALT2 and SiFTO lead to similar cosmological results. However, using the same 
dataset including SDSS-II, Ref. [90] discovered a difference of 0.2 in the best-fits of w by SALT2 and 
MLCS2k2, which already exceeds the total error budget (statistical and systematic) contributed from 
other sources (also see e.g. Ref. [95], where the discrepancy reaches 3cr). Furthermore, it is found 
that irrespective of which fitter is selected excluding MLCS2k2, different SNIa samples give largely 
consistent results, however as long as MLCS2k2 is considered, different datasets lead to notably dif- 
ferent constraints on Q.^^ [96]. The major difference between the MLCS2k2 fitter and the other two 
originates from the difference in the rest-frame U band region, where the training of MLCS2k2 relys 
exclusively on nearby SNIa observations while the other two also consider high redshift data [87] . 
In order to focus on the constraints on neutrino properties, we decide to use the largest sample inso- 
far (i.e. Union2. 1) which is obtained from the currently favored and more widely used method (i.e. 
SALT2). This also allows comparisons with most of other published results. 



'We directly take the best-fit values of a, b, and c presented in Table 2 of Ref. [82]. We have numerically checked 
that the correlation coefficients between these nuisance parameters and cosmological parameters are trivial enough to be 
neglected. The similar approach is adopted in Ref. [44]. 
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2.3 Baryon Acoustic Oscillations 



As the CMB radiation decouples from the primordial photon-baryon plasma, the features of acoustic 
oscillations are imprinted onto matter clustering as well and appear as peaks in the galaxy correlation 
function with a characteristic comoving separation of 100 h/Mpc [97, 98]. Baryon acoustic oscilla- 
tions (BAO) thus supply us with a standard ruler to measure distances out to the redshift where the 
bulk of galaxies are observed. 

The first detection was achieved by Ref. [99], using SDSS DR3 LRG sample with effective 
redshift z - 0.35. They defined the acoustic parameter which is independent of dark energy models 
as. 



Dwiz) 



ZC 



(2.12) 



where the distance combination is depicted by 
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with r{z) given by Eq. 2.1. We use the BAO measurements from the surveys of SDSS DRV [100] and 
WiggleZ [101]. Hence in total BAO results at five redshifts, i.e. z ^ 0.2, 0.35, 0.44, 0.6, 0.73, are 
involved in the constraints. Using the earlier SDSS result, Ref. [23] suggested that the inclusion of 
BAO can help break the degeneracy between dark energy and neutrino masses. 



2.4 Observational Hubble Parameter Data 

The differential-age technique directly measures the Hubble parameter as 

1 dz 

H{z) = -Y—i-, (2.14) 

and therefore endowed the OHD with special functionality as standard clocks. The quantity dz/df is 
usually determined by measuring the age differences between passively evolving galaxies with nearly 
the same spectroscopic redshift [50]. We use the most up-to-date datasets summarized by [51], which 
comprises nineteen OHD measurements over the redshift interval from 0.09 <z< 1.75. 



2.5 Further Data and Priors 

Complementary to late-time large scale probes (including WL) and other standard geometrical indi- 
cators (e.g. SNIa, BAO, OHD), the observation of CMB anisotropics is vital in this analysis. We 
use four independent power spectra, i.e. the temperature auto-power Cj^ , the curl-free component 
polarization auto-power Cf^, the curl component polarization auto-power Cf^, and the T-E cross- 
correlation Cj^, at the last scattering surface during the decoupling epoch (z ~ 1090). We use the 
WMAP seven-year data via including the standard pipeline for computing the likelihood, supplied and 
maintained by the WMAP collaboration ^. Apart from that, a top-hat prior of 10 Gyr < Iq <2Q Gyr 
on the cosmic age is employed. Last but not least, in all data combinations, we always impose the 
HST prior on the Hubble constant as Ho ^ 74.2 ± 3.6 km/s/Mpc [102]. 

^available at the LAMBDA website, http : //lambda . gsfc . nasa . gov/ 
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3 Constraints on Cosmological Parameters 



The cosmological parameters we use are presented in Table 1 . Our most generalized parameter space 
is composed of the "Vanilla+Extended" sets of parameters, represented by the following vector, 

P = {nkh^,nj)Mh\eA,T,n„]n{lO^^A,),{^my,N,«,w}). (3.1) 

In addition, the pivot value of the primordial power spectrum is taken to be - 0.05 Mpc"^ For 
simplicity, we assume a flat geometry and purely adiabatic initial conditions. We note that a non- 
zero curvature might loosen the constraint on ^ niy [103], but the WMAP seven-year result strongly 
confines -0.00133 < < 0.0084 (95% CL), which justifies our assumption. Moreover, we do not 
consider a running spectral index or any tensor contribution, for the constraint on neutrino masses is 
not relaxed even when they are allowed [104]. We also do not consider an evolving equation of state. 

Our code is a modified version of the publicly available package CosmoMC [105] ^. The global 
fitting to parameter vector P is achieved through the exploration of the corresponding multidimen- 
sional parameter space with the Markov Chain Monte Carlo (MCMC) technique. In realizing the 
MCMC approach, the Metropolis-Hastings algorithm is implemented to generate sets of chains con- 
taining sample points distributed according to the overall likelihood in the parameter space with 
top-hat prior probability distribution on each input parameter (see Table 1). 

Aiming at joint analyses of cosmological probes, we select five combinations from the afore- 
mentioned datasets (see Sect. 2), i.e. CMB, CMB-fWL, CMB-i-BAO-fOHD, CMB-i-BAO-HSNIa, 
and CMB-i-WL-i-BAO-i-OHD-i-SNIa, the last of which is often quoted as the full combination here- 
after. For all data combinations, the HST prior is always used. For a progressive and intensive 
exploration on the extended parameter set, there exist six scenarios for parameter combinations, i.e. 
Vanilla-i-2^ niy, Vanilla-i-A^eff , Vanilla-i-X mv+NeS, Vanilla-i-w-i-^ niy, Vanilla-i-w-i-A'^eff, Vanilla-i-w-i-2^ niy+NeS. 
As a consequence, we have attempted more than thirty runs in total, so as to sample the likelihood 
distribution in multivariate space for different combinations of observational data. For each run, eight 
chains are simultaneously generated using parallel computation, and they stop as soon as the criterion 
of Gelman and Rubin test (R statistics) is satisfied. Generally speaking, after convergence, our chains 
contain 10^ points each and guarantee R - I < 0.01. Then these chains are thinned and joined, left 
with more than 20000 points for the final constraints on each scenario for each data combination. The 
details of our results are reported as follows. 

3.1 Vanilla+i; niy 

We first explore the parameter space of Vanilla-F^ niy. Our main constraint results are quantitively 
summarized in Table 2. Fig. 1 displays the one-dimensional marginalized posterior distribution on 
some important parameters, while Fig. 2 shows the two-dimensional confidence contours revealing 
some crucial parameter degeneracies. Here as usual, we assume total degeneracy for the three basic 
neutrino eigenstates. The effective number of neutrinos is fixed as its standard value, i.e. A^gff - 3.046. 
The constraints are drawn under the context of the popular cosmological constant model (w = -1). 
Like many other papers on constraining ^ niy while fixing N^fi' and w, we found a 95% CL upper 
bound on the sum of neutrino masses, J] my < 0.476 eV, using the full combination. (Hereafter 
we neglect the unit of "eV" for brevity.) This is a 10% improvement over the formerly reported 
upper limit given by Ref. [45] and Ref. [46] using similar data. The WL evidently contributes to the 
overall constraints, pushing forward from ^ niy < 0.524 (CMB) to ^ My < 0.496 (CMB-i-WL), even 
better than the result given by CMB-i-BAO-i-SNIa. It reveals one of the greatest advantages of the WL 

'http : // cosmologist . info/cosmomc 
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technique that it measures the matter fluctuation amplitude very accurately, as seen from Table 2; erg is 
pinned down to the narrowest 68% confidence range for CMB+WL. Moreover, the OHD give similar, 
even slightly better results than the SNIa do niy < 0.486 vs. ^ /Wy < 0.518), considering the vast 
difference between the statistical size of the two data sample (19 vs. 580). This may be persuasive 
in designing next generation spectroscopic surveys in pursuit of large sample OHD measurements 
[106, 107]. From Fig. 2, we see a mild proportionality between ^ ruy and 0^ while a slight inverse 
proportionality between ^ niy and Hq is also seen. The inverse proportionality between ^ ruy and erg 
should be counted as the strongest correlation. Those degeneracies naturally explain why our best- 
fit values for Dm, Hq and erg are mildly larger, a little smaller, and unambiguously smaller than the 
WMAP seven-year recommended values. Except for the three derived parameters, constraints on the 
Vanilla parameter set retrieve nearly the same results as reported by the WMAP seven-year analysis. 
In particular, including additional probes does not result in much better constraints on n^, and the 
degeneracy between ^ ruy and n^, is marginal, since the variation of ?is changes the global shape of 
matter power spectrum whereas ^ my only damps P{k) at small scales. 

3.2 Vanilla+iVeff 

Then we test the scenario of Vanilla-i-A'^eff against each cosmological data combination. Likewise, the 
main results are presented in Table 3, Fig. 3 (one-dimensional) and Fig. 4 (two-dimensional). Here 
the global fitting is completed under the assumption of massless neutrino with ACDM cosmology. 
Using the full combination of cosmological probes, we obtain A/^eff - 3.271 ± 0.367, which rules 
out additional light degrees of freedom at nearly 95% CL. This result is in support of the standard 
prediction of big bang nucleosynthesis, and is primarily due to the contribution of the WL data (A/^eff = 
3.334 ± 0.496 for CMB-fWL, whereas A^eff = 4.340 ± 0.817 for CMB alone). The OHD also tend to 
favor the standard value. Including OHD, we obtain A/^eff = 3.722 ± 0.418, which is fully consistent 
with Ref. [51] (A^eff = 3.7 ± 0.4). Once again OHD give rather tighter bounds on the effective number 
of neutrinos than the SNIa do, and this time the difference between their constraining powers is 
more clearly revealed. The SNIa data seem to prefer the largest A'^eff best-fit value with the largest 
uncertainty. In this case, the inclusion of the SNIa data does not improve the constraint on either Hq 
or Qdm^^- Fig- 4 displays several parameter correlations. As reported by Ref. [61], we see the strong 
anticorrelation between A^eff and Qdm/j^- The degeneracy between A^eff and Hq is also remarkably 
strong. The propotional trend between A^eff and Q.^ is mild and disappears when the full combination 
is employed. These correlations arise from the role A^eff plays in the radiation content, which decides 
the epoch of matter-radiation equality as 

1 , ^mh^ ;2\-l /o o\ 

1 + Zea = = (Llyh ) . (3.2) 

^ Qr 1 +0.227 lA^eff ^ 

CMB measures precisely the present-day photon energy density to be D.yh^ = 2.469 x 10"^ for 
Tyo = 2.725 K [9]. Given the rather narrow range of Zeq = 3141;';|j^ [70], Eq. 3.2 thereby explains 
the correlations observed between A^eff> ^m, and Hq. The degeneracy between A'^Eff and Clrtuh^ is also 
extracted from Eq. 3.2, since Q.mh^ = ^bh^ + ^Duh^ and meanwhile the physical baryon density 
Q],h^ is tightly constrained by the CMB observations. Unlike its inverse proportionality with X "Jy, 
CTg shows clear proportionality with A'^eff ■ Because the observation of WL provides an accurate and 
independent measurement on erg, we thus expect that WL plays the most influential role in current 
constraints. The inclusion of SNIa is inadequate to deal with parameter degeneracies associated with 
A'^eff, and that is the reason for some of the poor constraints given by CMB-i-BAO-i-SNIa. Like the 
scenario of Vanilla4-^ niy, the constraints on the Vanilla parameter set are similar to the WMAP 
seven-year results, and the degeneracy between A'^eff and is broken when all data are included. 
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3.3 Vanilla+D my+Ne« 



We next investigate the case when both ^ my and A^gflf are set free, under the assumption of cosmolog- 
ical constant. Our main constraints are presented in Table 4, Figs. 5 and 6. Undoubtedly, the inclusion 
of Neff as a free parameter tremendously weakens the constraints on ^ my. Interestingly, when the full 
combination is considered, a robust 68% lower limit is revealed and our result is ^ my = OAll'^^ ^^^ 
at 68% CL. Nonetheless, considering that the best-fit of our current result is almost excluded at 95% 
CL by our former constraint when A^gff is fixed my < 0.476), this may illustrate our worry on the 
robustness of all stringent upper bounds on X 'Wy ever reported in previous literature. The constraints 
on Nsff are also relaxed due to the freedom on ^ my. Unlike former results, the extra neutrino species 
is slightly favored, as our best constraint reads N^fi- = 3.740 ± 0.446. The WL data still shows great 
promise as it refines the constraints from (2^ my < 1.515, A'^Eff - 5.729 ± 1.274, CMB alone) to 
(Yjmy < 1.393, A^eff = 4.308 ± 0.924, CMB+WL), mainly due to its precise constraint on erg. The 
combination of CMB+BAO-i-OHD(+HST) (Fig. 6) exerts more dominant influence on the parameter 
constraints, as the cyan and magenta contours are more alike. Besides the enormous information from 
the WMAP seven-year data, the OHD-i-HST delivers direct knowledge on the expansion history. In 
combination with the BAO distance indicator, this combination therefore is capable of successfully 
breaking several hard parameter degeneracies, e.g. my vs. Hq), (A^gff vs. Qm)- However even 
against the full combination, some correlations of parameters still remain. For instance, a larger ^ my 
value shows a preference for a smaller erg yet a slightly larger Q.^, whereas a larger A^eff favors a 
larger Qdm^^ and/or Hq. Most importantly, the mild proportionality between (J] my, N^s) reasonably 
points to the fact that a larger total neutrino mass is preferred by a larger effective number. As an 
increasing number of observations infer extra light degrees of freedom, varying ^ my and Nss si- 
multaneously will be a non-trivial consideration to derive constraints on both ^ my and A'^eff . As in 
previous scenarios, the constraints on the Vanilla parameter set are mostly unaffected. 



3.4 Vanilla+H'+2 my 

The previous stage of our analysis is constructed under the hypothesis of the cosmological constant 
being the underlying dark energy model. In what follows, we shall assume the dark energy component 
has a variable EOS, w. The other free extended parameter is ^ my. The main results for this scenario 
are shown in Table 5, Figs. 7 and 8. Owing to the inclusion of w as a free parameter, the best 
constraint on the sum of neutrino masses relaxes to ^ my < 0.627 (with no lower limit observed). 
In this scenario, the primary predominance the SNIa data holds over the OHD is demonstrated, i.e. 
the special importance of the SNIa data in constraining w. In practice, the SNIa impose much tighter 
constraint (w ^ - 1 .074 + 0.088) than the OHD (w = - 1 .240 ± 0. 182). Due to the degeneracy between 
w and ^ my, CMB-^BAO-i-SNIa thus presents better result my < 0.688) than CMB-i-BAO-i-OHD 
gives (Yj my < 0.819). Moreover, the combination of CMB-i-BAO-i-SNIa furnish an overally strigent 
constraint, as the green contours almost overlap with the cyan ones shown in Fig. 8. Besides, adding 
the present WL data does not help much at improving the constraints when w is allowed to vary freely. 
Although the technique of WL has been forecasted to be a promising tool to discriminate between 
dark energy and neutrino masses [41], currently speaking the WL data do not provide this discerning 
ability, albeit its constraint on erg is still incomparable. The anticorrelation between ^ my and erg is 
strong as usual, while and w are found to be highly correlated. Nevertheless, the degeneracies 
between ^ my, 0^ and Hq are mild. Especially, the inclusion of the SNIa data manages to break the 
degeneracy of {Q.^ vs. w). 
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3.5 Vanilla+H'+iVefr 



The next parameter scenario to explore is Vanilla+w+A'^eff. Table 6, Figs. 9 and 10 present the main 
results. UnUke the great impact varying w exerts on the ^ my constraints, the best constraint on A^gff is 
only compromised a little by freeing w, which reads A^eff - 3.454±0.386. But the inclusion of the free 
EOS parameter does make the combination of CMB(+HST) lose their capability to derive a reliable 
upper limit on NeS- However the importance of WL reappears unlike the case of Vanilla+w+X /My; 
the WL data constrains A^gff = 3.192 ± 1.214. Here the OHD wields incomparable constraining 
power on the effective number of neutrino species, i.e. A^eff - 3.623 ± 0.432, whereas the SNIa 
luminosity distance measurements still overpredict A^gff to an exceedingly large value. On w, the 
SNIa data without the contributions of WL or OHD strongly constrain w = -0.986 ± 0.077. In 
Fig. 10, interesting features again show up. The proportionality between Nes and Odm^^ is always 
strong and impossible to suppress. The degeneracy between A^eif and Q.^ is broken when the full 
combination is employed. The correlation between A^eff and Hq is present, yet weak. Particularly, 
the degeneracy directions between A^gflf and w, given by CMB+BAO+OHD and CMB+BAO+SNIa, 
show a nearly orthogonal feature. It is thereby enlightening since the OHD have the advantage of 
better constraining A^gff while the SNIa are capable of pinning down w more precisely. So when the 
OHD and the SNIa are combined, it is easy to break the degeneracy between A'^eff and w, as shown in 
our plots. 



3.6 Vanilla+H'+i; niy+Neff 

Finally, the most generalized scenario, where all of the three extended parameters are free to vary, 
is analyzed. Our main results are given in Table 7, Figs. 11 and 12. The best-fits give = 
0.556+°-23i, N^g = 3.839 ± 0.452, and w = -1.058 ± 0.088. Once again a 68% CL lower limit for 
X My is robustly observed. As discussed in previous scenarios, the WL data contribute most to the 
constraint on erg, and the narrowest 68% confidence intervals of w and Qm result from the inclusion 
of the SNIa data. In addition, due to the degeneracy of (X my vs. w), the SNIa offer better constraint 
on "^v than the OHD. However, concerning A'^eff , it is still the OHD that provide the most strict 
constraint. Moreover, it is found that for various scenarios, when A^gff is kept free, CMB+BAO+SNIa 
presents much weaker constraints on Quuh^ than CMB+BAO+OHD, yet their constraints look sim- 
ilar otherwise. This is a straightforward consequence of Eq. 3.2. Since CMB+BAO+OHD delivers 
direct information on the expansion history, it provides the most stringent constraint on Qdm^^- As 
seen from Fig. 12, when the entire set of extended parameters are freed, some of the parameter cor- 
relations turn mild, e.g. (A'^eff vs. Ho) and (Yj my vs. w). Nevertheless, the degeneracies of (2 niy vs. 
CTg) and (A^eff vs. Qdm^J^) are persistent. Once again the confidence contours of (Neff vs. w) display 
remarkable features: CMB+BAO+OHD and CMB+BAO+SNIa impose strict constraints with nearly 
orthogonal degeneracy directions, which emphasizes the importance of combining OHD and SNIa in 
future studies. Albeit mild, the proportion between ^ my and A^efif still exists. It again implies that 
one should vary ^ my and A^gff simultaneously in order to obtain more reliable constraints on cosmic 
relic neutrino properties. 

We note that within most scenarios, the best-fit values of erg are slightly low compared with the 
observations of galaxy clusters or Lya forest [21, 108]. However we argue that the results are strongly 
affected by the anticorrelation between X niy and erg. In our trial run with all the extended parameters 
fixed at their standard values, i.e. the usual ACDM scenario, we obtain erg = 0.791 ± 0.018 and 
fim - 0.255 ± 0.015, which are in total agreement with Ref. [83]. 
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4 Neutrino Impact on Structure Formation 



Over thirty sets of constraint results have been presented and discussed. In order to understand the 
physical origins underneath those parameter degeneracies, we investigate specifically the influences 
of the three extended parameters on growth of perturbation and formation of structures. We calculate 
a variety of matter power spectra via the CAMB code [109] ^ with the HALOFIT-y algorithm [80]. 
Here the default parameter set ("ALL") refers to the best-fit values of the Vanilla+w+^ my+Nes 
scenario against the full data combination (see Table 7). Hereafter, for simplicity, we assume total 
degeneracy for the standard three eigenstates of neutrinos, with mass for each flavor given by itiy = 

As mentioned in Sect. 1, massive neutrinos influence structure formation due to their large 
thermal velocity dispersion [110] 



15^(5) r,(z) 

\ ^(3) My 

where ^(s) denotes the Riemann zeta function, and Ty represents the neutrino temperature, which can 
be associated with the present-day temperature for the CMB, i.e. Ty{z) - (4/11)'^^ Tyo{l + z)- This 
kinematic activity, i.e. the efi'ect of free-streaming, prevents neutrinos from clustering below the free- 
streaming scale, and thus ensures the density perturbation of neutrinos below this scale negligible. 
By analogy with the Jeans length, this characteristic redshift-dependent scale is defined as [111], 

^Fs(z)- a/oTTT^^' (4.2) 
2(1 -i-z) Cs.y(z) 

with H(z) given by Eq. 2.3, and Cs^y being the sound speed for thermalized neutrinos, which can 
be strictly related to the velocity dispersion through Cs_v = ~y^v, in the non-relativistic limit [112]. 
Neutrinos become non-relativistic when their energy density (E) = J^^^yjTy = my, and this transition 
occurs at 1 -I- znr ~ 1890(my/l eV). Therefore for masses smaller than ~ 0.6 eV (the HALOFIT- 
y algorithm is well calibrated within this mass range), neutrinos will not dramaticaly slow down 
until the last scattering of photons ceases. Apart from that, if the neutrinos are massive enough, i.e. 
My > 0.005 eV, to have decelerated before dark energy dominates the expansion history, the Hubble 
horizon at the epoch when they turn into non-relativistic can be described by [1 10], 

knr - 0.0145 



1 +ZNR 



(^) </^h/Mpc. (4.3) 



This is the largest scale where the presence of finite-mass neutrinos starts to influence low-reshift 
matter power spectra by the free- streaming effect; above this scale, massive neutrinos can be treated 
simply as aggregating CDM. At scales k > /:nr> neutrinos not only free-stream to erase their fluctua- 
tions, but also damp the amplitude of the total matter power spectra by at least a few percent. 

We first display the signatures of massive neutrinos on matter power spectra at redshift z = 
0, 1,3, in Fig. 13. From Panel (a), as neutrinos become heavier, the suppression is increasingly 
severe at small scales, with the maximum suppression in the non-linear regime represented roughly by 

AP/P 9.8 (Qy/Qm) [113]. The redshift-dependence of the maximum suppression is analytically 

captured by the HALOFTT-y mapping algorithm [80]. From a crude glance at Panel (b), as more 
extra light degrees of freedom are introduced, matter power at small scales are more damped as well. 

"•http : // camb . info/ 
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Under more cautious inspection, the acoustic peaks are meanwhile shifted to smaller scales. We 
attempt to interpret these results in what follows. 

With a fixed total matter density Q.mh^, the increase of Nsfi' will postpone matter-radiation equal- 
ity (see Eq. 3.2). Since at subhorizon scales, the matter perturbations grow more efficiently during the 
matter-dominated epoch rather than during the radiation-dominated epoch, it is reasonable to deduce 
that the matter power spectrum is more severely damped at small scales relatively to lai^ge scales, 
which is one of the most pronounced features displayed in Fig. 13(b). For CMB, this effect is oppo- 
site; the height of the first acoustic peak is increased, mainly through the early integrated Sachs-Wolfe 
effect [1 14]. On the contrary, for the third and higher order peaks, the amplitudes are more damped 
due to the free-streaming of ultrarelativistic neutrinos [115]. The position of the «th acoustic peak 
can be estimated by /„ ~ nnjOA, with the acoustic scale parameter defined as (also see Table 1), 

- (4.4) 

KZdec) 

where r is given by Eq. 2.1 and Zdec ~ 1090 signals the epoch of photon-baryon decoupling [116]. 
Furthermore, Zdec is insensitive to the presence of massive neutrinos and r(zdec) almost remains unal- 
tered for different values of N^a [117]. Usually, the sound horizon is calculated via [73], 

rm^z) da 

Jo a^H{a) 

with the sound speed for the photon-baryon fluid c]. = c^/ [3(1 + R)] and R - (3Qb/40y)a. Since 
the massive neutrinos we are interested in maintain their relativistic speed before the recombination 
epoch, it is reasonable to conjecture that every neutrino species can be counted as a part of the radi- 
ation during that period. By increasing (decreasing) the value of A^eff > equivalently H{a) is changed 
according to Eq. 2.3, and the sound horizon at the time of decoupling, which dictates the position of 
the acoustic peaks, is thereby shrinked (enlarged). As a consequence, the acoustic peaks are shifted to 
smaller (larger) scales, corresponding to larger (smaller) wavenumbers of the matter power spectrum, 
particularly shown in Fig. 13(b). This feature, which is unique to the variation of A'^eff , might provide 
a smoking gun signature for the existence of extra light degrees of freedom. 

Unlike the impacts from (X my, Nsg), the effect introduced by varying w is a global change of 
magnitude for the matter power spectrum shown in Fig. 13(c). As inferred from Eq. 2.3, a larger 
value for w enables the dark energy to present itself with a non-vanishing amount in earlier stage of 
the expansion history. It thereby explains why the dotted curves in Fig. 13(c), corresponding to matter 
power spectra at redshift z = 3, are almost immune against the decrease of w deviating from -1. The 
change of w directly modifies the expansion history, it is therefore reasonable that the matter power 
at every scale is simultaneously enhanced/reduced. In addition, as shown, the redshift-dependence of 
maximum suppresion is still captured by the fitting algorithm of HALOFIT-y. Although the original 
HALOFIT mapping is only calibrated for the ACDM cosmological model, through appropriate mod- 
ifications of the growth function [1 18] or a specific fitting prescription developed by Ref. [1 19], one 
can reliably quantify the non-linear effect of varying w, accurate to a few percent. 

Fig. 14 presents the WL convergence power spectra as calculated from Eq. 2.8, with each panel 
corresponding to the panels of Fig. 13. No wonder the convergence power spectrum is also reduced 
by the increase of ^ niv and A^eff , or the decrease of w, since it is given by the integraion of matter 
power spectra along the line of sight. After integration, we have lost the unique signature for A'^eff- 
Fortunately, the impact from w still holds the different behavior of a global change, which can be 
used to differentiate from the effects introduced by ^ my and A^eff ■ The maximum suppression on the 
convergence power is at / ~ 10^, insensitive to the variations of the three parameters. Besides, we 
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also plot the result predicted by the linear perturbation theory in Fig. 2.8(a). The primary discrepancy 
takes place below the scale corresponding to / > 10^, and it thereby illustrates the importance of an 
accurate knowledge of the full non-linear matter power spectrum, for future WL surveys to achieve 
their full potential [120]. 

Eventually, the numerical results of the aperture-mass variance given by Eq. 2.9 are demon- 
strated in Fig. 15, together with the actual data of CFHTLS-T0003. From Panel (a), we see that 
the current detection of WL signals mainly lies within the semi- to non-linear regime, where the 
non-Unear modeling of matter clustering is vital for a precise interpretation of the observational data. 
As repeatedly mentioned, the HALOFIT-v algorithm has great advantages itself, yet the effects of 
baryons are not covered. The hot baryons harbored within intracluster medium have been shown to 
affect the WL shear power spectrum up to ten percent at k ~ 10 h/Mpc [121]. If sufficiently strong 
AGN feedback is included, ten percent is already reached at ^ ~ 1 h/Mpc [122]. Using hydrodynamic 
simulations with consideration of baryon physics, Ref. [123] discovered that as much as 40 percent 
bias in the constraint on dark energy EOS can be introduced by ignoring the effects of baryons. How- 
ever, most of the baryon contribution is at smaller scales than currently observed. Besides, finite-mass 
neutrinos exert quite different impacts on matter clustering compared with the baryons. Thus using 
current optimal mapping algorithm, HALOFIT-v, will suffice for the constraints on ^ niy. In the 
precision analysis of the WL, not only are the theoretical prescriptions crucial, but the data quality 
is also influential. As shown in Fig. 15, a small bump is present around the semi-nonlineai" scales 
in the CFHTLS data. Although (^ap) ^ more localized filter and thus less prone to large scale 
systematics or E-/B-mode mixing, this feature is still mainly caused by some thorny systematic is- 
sues, e.g. shape measurements, intrinsic alignments, PSF anisotropics, catastrophic photo-z error. 
For a thorough discussion on WL systematics, observables and their interplay, we refer the reader to 
Ref. [124]. Anyhow, the community is looking forward to better settlement of these systematics, in 
the final release of CFHTLS lensing data in the near future. 

5 Conclusion and Discussion 

In this paper, we present a comprehensive analysis of the constraints on the sum of neutrino masses, 
the effective number of neutrino species, and the constant EOS parameter of dark energy. Via em- 
ploying the most recent observational data of CMB-FWL-i-BAO-i-OHD-i-SNIa, we have pinned down 
the total neutrino mass to ^ /My < 0.476 at 95% CL, whereas the constraint reads ^ my = OAll'^^ ^^^^ 
(68% CL) when N^^s is simultaneously set loose, and even gives ^ niy = 0.556;|;Q 2gg (68% CL) if w is 
also included as a free parameter. Given that we only consider the constant EOS for the dark energy 
in present work, the results including evolving dark energy component will provide even less strin- 
gent constraints. This highlights our concern that in constraining neutrino masses using cosmological 
probes, one needs to take into consideration the correlations between (^ my, Ns«, w) seriously. For 
the effective number parameter, we obtained the best constraints as A^gff = 3.217 ± 0.367 (with fixed 
X my and w), A^eff = 3.740 ± 0.446 (with freed niy and fixed w), A^eff = 3.454 ± 0.386 (with fixed 
Yj my and freed w), and A'^eff - 3.839 ± 0.452 (with both ^ niy and w freed). For the most general- 
ized parameter scenario, the standard model of N^^s = 3.046 is compatible at 1.75o" CL. If statistical 
uncertainties were dominating, which is probably true for the BAO and OHD, then we can rule out 
the standard three flavor neutrinos model at 5cr CL, when the future datasets over eight times as large 
as current samples are obtained. But for some probes like SNIa, systematics are already dominat- 
ing, so the expected improvements are not easily estimated. Most probably in the future, equivalent 
or improved constraints will be obtained with BAO angular diameter distances and Hubble parame- 
ter measurements. The constraints on the EOS parameter all along favor the cosmological constant 
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model, regardless of whether or not varying ^ niy and/or N^s. That is mainly due to the effectiveness 
of the SNIa data. Moreover, the OHD still provide helpful improvements on dark energy constraints. 
We also pay particular attention to parameter degeneracies. The strong correlation between A^eff and 
Hq (also A^eff and Qdm'^^) can be explained by the rather certain epoch for matter-radiation equality. 
The strong anticorrelation found in my vs. erg) results from the suppression of matter perturbation 
growth ascribed to the presence of finite-mass cosmic neutrinos. As shown in Sect. 3, due to the 
different proportionality trends that ^ niy and N^^s hold with respect to erg, as well as the advantage of 
WL in measuring erg, the technique of WL does show great potential in discriminating between the 
effects of Yj and A^eff for future large surveys. 

Another key aspect in our work is the investigation on the constraining capabilities of different 
cosmological probes. The WL observation at current stage does contribute to some improvement 
if the dark energy parameter is fixed. However when w is free to vary, adding the WL data does 
not result in much better constraints on ^ my. Perhaps the WL tomography and its cross-correlation 
with other tomographic probes will help improve this situation [42, 43]. On the comparison between 
the results from the inclusions of OHD and SNIa, it is found that for the constraints on Yj ^v, when 
dark energy EOS is fixed (w - -1), the result given by OHD is similar to that by SNIa luminous 
distances or BAO angular diameter distances, with OHD slightly better. If N^^s is set free, OHD 
distinctly prevail over SNIa. Nevertheless, as long as w is freed, the SNIa data confine Yj "^v more 
strictly than the OHD do. Concerning the constraints on the Hubble constant, when only Y is 
set loose, the SNIa data provide similar constraints as the OHD do. Once N^s is varying, the OHD 
perform much better than the SNIa, due to the overwhelming power of the OHD in constraining N^s- 
Additionally, the SNIa and the OHD provide almost the same constraints on Q.^. On the degeneracy 
between N^s and w, the inclusions of OHD and SNIa lead to nearly orthogonal degeneracy directions, 
as explained by their specific powers in constraining A^gflf and w, respectively. As a consequence, their 
union will facilitate the breaking of this parameter degeneracy so well that our global constraints on 
all cosmological parameters will benefit. 

The tantalizing indication of a non-zero total neutrino mass is really interesting (see Sec. 3.3 
and Sec. 3.6). We argue this is the result of the combined effect from all cosmological probes in 
our study. First of all, due to the large parameter space and the fact that some parameters, e.g. 
Qb/i^ and r, are badly constrained by other probes, the CMB data are indispensable, yet not enough 
even when combined with the Hq prior. The WL technique provides plenty of information on mass 
distribution and matter clustering, which are sensitive to the neutrino total mass, however its use 
is limited by measurement errors as well as the inclusion of variable w. Although not directly re- 
lated to neutrino masses, the late-time probes, i.e. BAO, SNIa and OHD, are particularly helpful 
in disentangling some parameter degeneracies. Thereby combining all the observations together, 
i.e. CMB-i-WL-i-BAO-i-OHD-i-SNIa(-i-HST), we can constrain Y "Jy very precisely and even impose a 
68% CL lower limit. As also pointed out recently by Ref. [125], neither future CMB experiment or 
weak lensing shear measurement or galaxy tomographic survey alone can reach ~ 5cr detection of a 
non-zero sum of neutrino masses; it is their joint analysis that effectively breaks parameter degenera- 
cies. 

The different impacts that dark energy, cosmic neutrino total mass and effective number exert 
onto matter clustering and structure formation are also discussed. Regarding the effects on low- 
redshift matter power spectra, the increase of w tends to reduce the global magnitude of matter power 
at all scales, due to direct modification of the expansion law. On the contrary, the change in Y 
and/or A^gff only affects the power spectrum at scales k > /cnr. The increase of Y "iv damps more 
severely the non-linear matter power due to the free-streaming effect, while the increment on A^eff 
causes delay to the matter-radiation equality epoch, and thus suppresses the growth of matter pertur- 
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bation at small scales. The suppresion caused by Ngg is nearly redshift-independent for < 1 h/Mpc. 

In addition, we note an unique feature of A^eff : its variation is capable of shifting acoustic peaks 
in late-time matter power spectra. Accordingly, the WL convergence power spectrum is altered, and 
the aperture-mass variance affected. It is also discussed that the precision analysis of WL demands 
the accurate knowledge of the full non-linear matter power spectrum and the next stage cosmological 
simulations had better take gas dynamics into account. Albeit the current data of CFHTLS still suffer 
from systematics, high hope and great endeavor have been invested in this program, which will finally 
promote our knowledge on the formation and evolution of cosmic structure through WL. 

Recently, Ref. [126] tested the model of eV mass scale sterile neutrinos, employing two differ- 
ent SNIa datasets given by different light curve fitters, i.e. Union2 with SALT2, and SDSS-II with 
MLCS2k2. They found that such sterile neutrino model is strongly supported by the SNIa sample of 
SDSS-II with MLCS2k2 (~ 55 times more probable than the null model), whereas disfavored by the 
Union2 SNIa data with SALT2, as well as other cosmological probes. This is not entirely unexpected 
since several studies have already shown that MLCS2k2 seems to support more exotic physics rather 
than the ACDM model, whereas the very same data standardized with SALT2 are fairly consistent 
with cosmological constant as well as other types of cosmological observations [96, 127, 128]. We 
note that in judging the viability of Ref. [126]'s sterile neutrino model, the different fitters do play 
a role, but also the different SNIa data samples and the estimates of the systematic errors associated 
with the samples. Therefore to conclude that the light curve fitter is the key factor of vindicating 
such sterile neutrino model or not, more work should be dedicated to the re-fitting of the same SNIa 
dataset, especially Union2.1, using different methods. 

Cosmological probes we have discussed are mainly sensitive to the sum of neutrino masses, 
yet almost blind to mass splitting scenarios, if any, for different neutrino eigenstates. Fortunately, 
cosmological observations can reliably rule out total degeneracy between the eigenstates if 2 f^v is 
constrained robustly to be less than 0.22 eV, and furthermore nail down the Normal Hierarchy scenario 
as the true mass splitting when ^ niy < 0.1 eV is conservatively derived [cf. Ref. 129, Figure 1]. 

However, as shown in Fig. 15, these small values for X"^v are seemingly disfavored by the 
M^p^ data alone, as the constraint on A^eff tends (slightly) to prefer additional light degrees of free- 
dom. In practice, there have been some arguments on discriminating the neutrino mass hierachy 
within the framework of Fisher formalism [130-133]. Ref. [134] conducted N-body simulations to 
study seperately the effects from different mass hierarchies and obtained a ~ 0.5% difference. This 
small, yet non-trivial result underlines its measurable potential from upcoming space surveys (e.g. 
Euclid), in combination with ground-based experiments. In our present analyses, A^eff is treated as the 
effective number of neutrinos. We do realize that except members of the neutrino family, other more 
exotic candidates, e.g. axions, gravitinos, or even primordial magnetic fields, can equally fill for this 
extra relativistic component. Besides, early dark energy has been shown to be able to reconcile the 
standard model of three families of neutrinos with various cosmological observations [67]. Nonethe- 
less, as large scientific programs are already providing new results, we are confident that the truth of 
dark radiation will soon be unveiled. 
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Figure 1. One-dimensional marginalized constraints on individual cosmological parameters from five 
combinations of different obsei-vations, i.e. CMB (blue), CMB+WL (red), CMB+BAO+OHD (magenta), 
CMB+BAO+SNla (green), and CMB+WL+BAO+OHD+SNla (cyan). Here 2 niy is free whereas A^eff and w 
are fixed to theii- standard values (A^eff = 3.046, w - -1). 



Table 1. Key parameters and their prior information. 



Set 


Parameter 


Symbol 


Prior 


Vanilla 


Physical baryon density 




[0.005,0.1] 




Physical dark matter density 




[0.01,0.99] 




100 X Angular size of sound horizon 


100 X 9a 


[0.5, 10] 




Reionization optical depth 


T 


[0.01,0.8] 




Scalar spectral index 




[0.5, 1.5] 




Scalar spectral amplitude 


In (10^%) 


[2.7,4] 


Extended 


Sum of neutrino masses 


Zm, [eV] 


[0, 10] 




Effective number of neutrinos 


A^eff 


[0, 10] 




Constant dark energy EOS 


w 


[-2,0] 


Derived 


Total matter density 








Power spectrum normalization 








Hubble constant 







The "Vanilla" set of parameters is the basic set that is always fed into our MCMC code with 
uniform prior distribution, as shown within the table. What we are exceedingly interested in 
is the "Extended" set of parameters, whose flat prior distributions are also given. By fixating 
some members of this set while varying the rest, it is possible to investigate individually and 
collectively the constraints they receive from observational data. And meanwhile we are able to 
learn some knowledge about how parameter degeneracies affect the constraints. As shown in this 
table, Qm, erg and Hq are always derived. Additionaly, a flat prior ([0,2]) is assigned to Agz, i.e. 
the SZ template amplitude, which is always marginalized as a nuisance parameter. 
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Figure 2. Two-dimensional confidence contours of several parameter pairs from five combinations of different 
obsei-vations, i.e. CMB (blue), CMB+WL (red), CMB+BAO+OHD (magenta), CMB+BAO+SNIa (gi'een), 
and CMB+WL -HBAO+OHD+SNIa (cyan). For each data combination, the inner and outer contours refer to 
the borders of 68% and 95% confidence regions, respectively. As in Fig. 1, 2 is free here, but A^gft and w 
are fixed. 
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Figure 3. Same as Fig. 1, except that A^eff is a free parameter here, with massless neutrinos and w = -1. 
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Figure 6. Same as Fig. 2, but for both 2 '«v and Neff kept free, with w - -I. For the sake of clarity, the 
dotted lines in each panel denote the same confidence regions for CMB+BAO+OHD, as shown by the magenta 
contours. 
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Figure 7. Same as Fig. 1, but also with w freed, and A^eff is still fixed. 



-29- 





-30- 




Figure 12. Same as Fig. 2, yet all members of the extended set (J] niy, Neg, w) are treated as free parameters. 
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Figure 13. Matter power spectra and their fractional differences. In each panel, the sets of solid, dashed 
and dotted lines represent power spectra at redshift z - 0, 1,3, respectively. The blue curves in each panel 
are calculated from the best-fit parameter values when all cosmological probes are considered, with all three 
extended parameters (2 '«v> ^eif> m') kept free (see Table 7). Compared with the blue set of lines, other sets 
are computed with only one given parameter altered to the values shown within each panel. 
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Figure 14. Convergence power spectra and their fractional differences. As in Fig. 13, the blue curve corre- 
sponds to the best-fit parameter values given all cosmological data with all of (2 niy, A^eff, w) free (see Table 7). 
Other lines are plotted with only one given parameter changing its value in order to show the effects of the three 
extended parameters. 
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Figure 15. Aperture-mass variance from observation (CFHTLS) and from theoretical calculations. As in 
Fig. 13, the blue curve is given by the best-fit parameter values against the full combination with all of 
(Yj my, Nss, w) free (see Table 7). Other lines are drawn with only one given parameter altered. 



Table 2. Constraints on Parameter Scenario of "Vanilla+2 niy' Using Multiple Combinations of Cosmological 
Probes. 







CMB 


CMB+WL 


CMB+BAO+OHD 


CMB+BAO+SNIa 


CMB+WL+BAO+OHD+SNIa 


Vanilla 


lOOfib/!- 
IOODdm/i^ 

T 

«j 

ln(10'%) 


2.242 ± 0.055 
11.04 ±0.49 
103.91 ±0.25 
0.090 ±0.015 
0.970 ±0.013 
3.073 ± 0.036 


2.241 ±0.051 
11.00 ±0.35 
103.90 ±0.25 
0.089 ±0.015 
0.969 ±0.013 
3.070 ±0.031 


2.242 ±0.051 
11.33 ±0.32 
103.94 ± 0.25 
0.088 ± 0.014 
0.968 ± 0.013 
3.081 ± 0.034 


2.236 ± 0.052 
11.33 ±0.31 
103.91 ±0.26 
0.087 ±0.015 
0.967 ±0.012 
3.076 ±0.035 


2.229 ± 0.050 
11.19±0.24 
103.88 ±0.23 
0.087 ±0.014 
0.964 ±0.011 
3.068 ±0.030 


Extended 


2m„ [eV] 


< 0.524 


< 0.496 


< 0.486 


< 0.518 


< 0.476 


Derived 


/fo/100 


0.272 ± 0.028 
0.755 ± 0.046 
0.700 ± 0.024 


0.270 ± 0.023 
0.751 ±0.035 
0.702 ± 0.022 


0.286 ±0.015 
0.766 ± 0.048 
0.689 ±0.012 


0.288 ±0.015 
0.758 ± 0.049 
0.687 ±0.013 


0.283 ±0.013 
0.743 ± 0.035 
0.689 ±0.012 



Concerning the reported results of all parameters except 2 "V^ the best-fit values refer to the means calculated from 
one-dimensional marginalized posterior probability distributions, as plotted in Fig. 1 for some parameters. The reported 
errors originate from the symmetric 68% confidence intervals with respect to the corresponding means. For 2 m,, the 
95% upper limit is given. Here N^ff and w are fixed to their standard values (A'eff = 3.046, w = -1). 



Table 3. Constraints on Parameter Scenario of "Vanilla-i-A^eff" Using Multiple Combinations of Cosmological 
Probes. 







CMB 


CMB+WL 


CMB+BAO+OHD 


CMB+BAO+SNIa 


CMB+WL+BAO+OHD+SNIa 


Vanilla 


lOO^DM^^ 

lO-^^A 

T 

In (10'%) 


2.215 ±0.057 
13.54 ± 1.85 
103.42 ± 0.38 
0.086 ±0.015 
0.978 ±0.014 
3.126 ±0.048 


2.228 ± 0.057 
11.15 ±0.76 
103.74 ±0.37 
0.089 ±0.015 
0.971 ±0.014 
3.074 ± 0.032 


2.214 + 0.054 
12.56 + 0.83 
103.57 ±0.32 
0.084 ± 0.014 
0.969 ± 0.012 
3.105 ±0.035 


2.213 + 0.052 
14.50 ± 1.50 
103.30 ±0.35 
0.084 ±0.014 
0.982 ±0.014 
3.139 ±0.040 


2.204 ± 0.052 
11.39 ±0.57 
103.64 ±0.33 
0.081 ±0.013 
0.961 ±0.011 
3.064 ± 0.029 


Extended 


A'eff 


4.340 ±0.817 


3.334 ± 0.496 


3.722 ± 0.418 


4.827 ±0.843 


3.271 ±0.367 


Derived 


o-g 
//o/lOO 


0.283 ±0.031 
0.870 ± 0.053 
0.746 ± 0.026 


0.250 ±0.013 
0.799 ± 0.023 
0.732 ± 0.024 


0.287 ± 0.015 
0.848 ± 0.029 
0.717 ±0.016 


0.289 ±0.015 
0.892 ±0.041 
0.760 ± 0.032 


0.270 ± 0.009 
0.804 ±0.018 
0.710 ±0.016 



Same as Table 2, except that N^ft is free to vary, with massless neutrinos and w = -1. We present the results of A'efr 
following the usual convention as introduced in Table 2. 
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Table 4. Constraints on Parameter Scenario of "Vanilla+2 niy+Ngg' Using Multiple Combinations of Cosmo- 
logical Probes. 







CMB 


CMB+WL 


CMB+BAO+OHD 


CMB+BAO+SNla 


CMB+WL+BAO+OHD+SNIa 


Vanilla 


lOOfib/!- 


2.186 ±0.058 


2.204 ± 0.059 


2.214 ± 0.055 


2.212 ±0.054 


2.213 ±0.055 




IOODdm/i^ 


16.75 ±2.95 


13.32 ± 1.89 


12.68 ±0.86 


14.35 ± 1.42 


12.34 ±0.80 






103.11 ±0.34 


103.43 ± 0.36 


103.56 ±0.32 


103.35 ±0.35 


103.58 ±0.32 




T 


0.084 ±0.014 


0.086 ±0.014 


0.088 ± 0.015 


0.087 ±0.014 


0.087 ±0.014 




«j 


0.985 ±0.014 


0.976 ±0.014 


0.974 ±0.013 


0.987 ±0.015 


0.972 ±0.012 




In (10'%) 


3.157 ±0.050 


3.097 ± 0.038 


3.101 ± 0.037 


3.132 ±0.041 


3.086 ±0.033 


Extended 


2m.-[eV] 


< 1.515 


< 1.393 


< 0.758 


< 0.775 


Q 421+0.186 




A'eif 


5.729 ± 1.274 


4.308 ± 0.924 


3.843 ± 0.439 


4.814 ±0.889 


3.740 ± 0.446 


Derived 




0.341 ±0.051 


0.297 ± 0.040 


0.296 ±0.018 


0.294 ±0.015 


0.292 ±0.016 






0.806 ± 0.063 


0.729 ±0.051 


0.762 ± 0.059 


0.807 ± 0.063 


0.731 ±0.037 




/fo/100 


0.746 ± 0.034 


0.724 ± 0.032 


0.709 ± 0.017 


0.749 ±0.031 


0.707 ±0.017 



Same as Table 2, but for both 2 my and N^s kept free, with w = —\. The conventions for presenting the results of 
2 niy and Nts are as usual. Only for the full data combination (CMB+WL+BAO+OHD+SNIa), the result of 2 riiy is 
represented by the best-fit value with asymmetric errors marking the ranges of 68% confidence interval. 



Table 5. Constraints on Parameter Scenario of "Vanilla+w+ J] niy" Using Multiple Combinations of Cosmo- 
logical Probes. 



Vanilla 



CMB 



CMB+WL 



CMB+BAO+OHD CMB+BAO+SNIa CMB+WL+BAO+OHD+SNIa 



lOO^DM^ 
lO^^A 
T 

In (10'%) 



2.191 + 0.059 
11.87 ±0.70 
103.79 ±0.26 
0.085 ± 0.014 
0.955 ± 0.015 
3.075 ± 0.034 



2.186 ± 0.062 
11.59 ±0.57 
103.74 ±0.26 
0.086 ±0.014 
0.954 ±0.015 
3.063 ± 0.030 



2.214 ±0.054 
11.82 ±0.45 
103.87 ±0.25 
0.085 ±0.014 
0.959 ±0.013 
3.081 ±0.034 



2.218 ±0.052 
11.53 ±0.43 
103.85 ± 0.26 
0.087 ±0.015 
0.962 ±0.013 
3.077 ± 0.036 



2.225 ± 
11.29 = 
103.88 
0.086 ± 
0.964 ± 
3.067 ± 



0.054 

:0.31 

tO.24 
0.014 
0.012 
0.031 



Extended Urn, [eV] < 1.079 < 1.072 < 0.819 < 0.688 

w -1.379 ± 0.247 -1.340 ± 0.237 -1.240 ±0.182 -1.074 ±0.088 



< 0.627 
- 1.034 ± 0.080 



Derived 



/fo/100 



0.255 : 
0.788 : 
0.746 : 



0.029 
0.068 
0.038 



0.250 ± 0.028 
0.757 ± 0.047 
0.746 ± 0.038 



0.268 ±0.018 
0.794 ± 0.065 
0.725 ± 0.029 



0.286 ±0.016 
0.762 ± 0.062 
0.693 ±0.017 



0.283 ±0.014 
0.741 ±0.035 
0.692 ±0.016 



Same as Table 2, but also with w freed, and N^fs is still fixed. The results of w is presented in the conventional manner. 
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Table 6. Constraints on Parameter Scenario of "Vanilla+w+A^efF" Using Multiple Combinations of Cosmolog- 
ical Probes. 



CMB CMB+WL CMB+BAO+OHD CMB+BAO+SNIa CMB+WL+BAO+OHD+SNIa 



Vanilla 


lOOQ.t,h- 


2.227 ± 0.061 


2.230 ± 


0.061 


2.207 ± 0.053 


2.215 ± 


0.054 


2.212 ± 


0.054 




IOODdm/?^ 


17.09 ± 3.72 


11.02 d 


I 1.61 


12.52 ±0.84 


14.05 d 


: 1.56 


11.50d 


:0.57 






103.16 ±0.43 


103.82: 


±0.49 


103.60 ± 0.32 


103.36: 


±0.35 


103.60: 


±0.32 




T 


0.086 ±0.015 


0.087 ± 


0.015 


0.084 ±0.014 


0.085 ± 


0.014 


0.084 ± 


0.014 






1.002 ± 0.028 


0.965 ± 


0.026 


0.967 ± 0.014 


0.980 ± 


0.016 


0.967 ± 


0.013 




ln(10'%) 


3.173 ±0.060 


3.062 ± 


0.047 


3.107 ±0.036 


3.133 ± 


0.042 


3.071 ± 


0.030 


Extended 


A'eff 


> 2.885 


3.192 ± 


1.214 


3.623 ± 0.432 


4.608 ± 


0.860 


3.454 ± 


0.386 




w 


-0.794 ± 0.240 


- 1.076 H 


t 0.233 


-1.054 ±0.130 


-0.986 = 


t 0.077 


-0.937 H 


t 0.060 


Derived 




0.354 ± 0.078 


0.245 ± 


0.044 


0.282 ±0.019 


0.291 ± 


0.016 


0.278 ± 


0.013 




0"8 


0.849 ± 0.066 


0.807 ± 


0.050 


0.865 ± 0.056 


0.876 ± 


0.048 


0.785 ± 


0.025 




/fo/100 


0.742 ± 0.037 


0.740 ± 


0.037 


0.723 ± 0.026 


0.748 ± 


0.033 


0.703 ± 


0.018 



Same as Table 2, but here A'eff and w are instead set free, with massless neutrinos. As illustrated in Fig. 9, the upper 
limit of A'eff for CMB is severely affected by the preempted prior (A'eff < 10), therefore not trustable. Thus we present 
its 95% lower limit instead. 



Table 7. Constraints on Parameter Scenario of "Vanilla+w+2 '«v+^etf" Using Multiple Combinations of Cos- 
mological Probes. 



Vanilla 



CMB 



CMB+WL 



CMB+BAO+OHD CMB+BAO+SNIa CMB+WL+BAO+OHD+SNIa 



~l0(XVr~~ 
IOODdm^^ 

lO^^A 



T 

n. 



2.191 ± 0.062 
16.40 ±3.61 
103.24 ± 0.45 
0.085 ± 0.015 
0.985 ± 0.028 



2.191 ± 0.061 
13.35 ± 2.40 
103.51 ±0.47 
0.084 ± 0.015 
0.967 ± 0.024 



ln(10"'As) 3.146 ± 0.062 3.088 ± 0.046 



2.197 ±0.057 
13.08 ± 0.93 
103.53 ±0.31 
0.086 ± 0.014 
0.966 ± 0.014 
3.104 ±0.036 



2.211 ±0.055 
14.07 ± 1.51 
103.37 ±0.36 
0.087 ±0.014 
0.980 ±0.015 
3.126 ±0.041 



2.202 ± 0.056 
12.70 ±0.86 
103.52 ±0.32 
0.086 ±0.014 
0.970 ±0.013 
3.089 ± 0.032 



0.556 g 

3.839 ±0.452 
-1.058 ±0.088 



Extended 2™. [eV] < 1.606 < 1.539 < 1.035 < 0.869 

A'eff > 2.642 4.125 ± 1.426 3.794 ± 0.447 4.574 ± 0.896 

w -1.051 ± 0.335 -1.155 ±0.278 -1.206 ±0.183 -1.038 ±0.088 



Derived 



//o/lOO 



0.338: 
0.779 : 
0.745 : 



0.075 
0.070 
0.038 



0.292 ± 0.057 
0.729 ± 0.058 
0.733 ± 0.038 



0.282 ±0.021 
0.790 ± 0.067 
0.737 ±0.031 



0.295 ±0.016 
0.802 ± 0.065 
0.743 ±0.031 



0.294 ±0.016 
0.723 ± 0.036 
0.712 ±0.020 



Same as Table 2, yet all members of the extended set (2 niy, A'eff, w) are treated as free parameters. Likewise, for CMB, 
the 95% lower limit is given for A'eff. Also the result of 2 for the full combination is represented by the best-fit along 
with 68% confidence region. Note that the case where the entire extended set is kept free against the full combination 
is quoted as "ALL" in Section 4. 
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